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Abstract 

We use the semi-analytical program RCFORGV to evaluate radiative corrections to one-photon 
radiative emission in the high-energy scattering of pions in the Coulomb field of a nucleus with 
atomic number Z. It is shown that radiative corrections can simulate a pion polarizability effect. The 
average effect is a™ = —f3^ c = (0.20 ±0.05) x 10 _43 cm 3 , for pion energies 40-600 GeV. We also study 
the range of applicability of the equivalent photon approximation in describing one-photon radiative 
emission. 



1 Introduction 



Pion electric a and magnetic j3 Compton polarizabilities [l]-[4] characterize the induced tran- 
sient dipole moments of pion subjected to external oscillating electric E and magnetic H fields. 
The polarizabilities can be obtained from precise measurements of the gamma-pion Compton 
scattering differential cross section. They probe the rigidity of the internal structure of the 
pion, since they are induced by the rearrangement of the pion constituents via action of the 
photon electromagnetic fields during scattering. For the 771 interaction at low energy, chiral 
symmetry provides a rigorous way to make predictions. This approach [5, 6] yields a t = -/3 n 
= 2.7 ± 0.4 (in units 10" 43 cm 3 ). 

The radiative scattering of high energy pions in the Coulomb field of a nucleus [7]: 

tt( Pi ) + Z{p) - 7V( P2 ) + Z{p') + 7 (fc'), (1) 

determines the 7 7r Compton scattering: 

7r ( Pl ) +7 ( A; )^ vr (p 2 ) +7 (A;'). (2) 

The 77T scattering was measured [8, 9] with 40 GeV pions at Serpukhov via the Primakoff 
radiative scattering process (1). The incident pion scatters from a virtual photon, characterized 
by 4- momentum k = p—p'; and the final state gamma ray and pion are detected in coincidence. 
The pion electric polarizability a n was deduced in a low statistics (~ 7 000 events) experiment 
to be: 

a w = -Ar = 6.8 ± lA stat ± 1.2 syst . (3) 

It was assumed in the analysis that a n + (3 W = 0, as expected theoretically [5]. This result 
differs from the chiral prediction by more than two standard deviations. The experimental 
situation points to the need for much higher quality data and more attention to the systematic 
uncertainties arising from different measurement and analysis techniques. A new Primakoff 
experiment is planned at FNAL [10]. 

In the lowest order, the process of the radiative scattering of the pions in the Coulomb 
field of a nucleus is described by Feynman graphs of Fig. 1: 1 la) the QED interaction of the 
pointlike pion with the electromagnetic field of the nucleus; Id) the two-photon interaction of 
the pion, including the structure effects described by the polarizabilities of the pion. In the 
next order of a, there are contributions from the QED processes presented in Fig. l(b,c). The 
double bremsstrahlung process (lc): 

tt(pi) + Z{p) - Trip 2 ) + Zip') + 7 (fci) + 7 (fc 2 ), (4) 

and the one- loop diagrams (lb) with the virtual photon, both contribute to the measured 
cross section of the reaction (1). 

The total cross section for the radiative scattering process (1) can be written in the form: 

da tot = da Born + do voi + da rc . (5) 

Here, da Born corresponds to scattering of a pointlike pion in the Born approximation, and 
the following terms give the contributions of the structure effects (linear terms in the pion 
polarizability) and the radiative corrections. The radiative corrections are important for high 
statistics experiments, and are needed to extract the Compton cross sections from data. The 



1 



a) 


7 


b) 




Z 4 


7 

! ► 














d) 


— ► 






w w 



Figure 1: Feynman graphs for radiative scattering of pions on Coulomb center: a) the Born 
approximation; b) the one-loop QED corrections; c) the double photon bremsstrahlung; ) the 
structure effects - pion polarizabilities. 

experimental data can be interpreted in terms of polarizabilities, only after accounting for the 
radiative corrections. 

The theoretical investigation [11]- [13] of these radiative corrections (RC) was undertaken 
initially for the analysis of the experiment [8, 9] of JINR (Dubna) and IHEP (Serpukhov). 
It was shown that the detailed properties of the gamma detector are important; such as 
the photon detector threshold, energy and position resolution, and the two-photon angular 
resolution. For the setup of Antipov et al. [8, 9] at 40 GeV, the radiative corrections were 
estimated to affect polarizability determinations at the level of ~ 3% [13]. 

In this note we present a study of the radiative corrections to the process (1) at higher 
energies and for higher Z targets. Our calculations were carried out with the semi-analytical 
program RCFORGV [14]. 

2 The semi-analytical program RCFORGV 

The analytical formulae for the RC to the inclusive spectra of the scattering of a spinless point- 
like particle by a Coulomb center with emission of a hard photon were obtained in [11, 13]. 
The contribution of the QED one-loop corrections (Fig. lb) and of the double bremsstrahlung 
from pions (Fig. lc) were calculated exactly without any approximation and for m n ^ 0. The 
calculation of the Feynman diagrams was carried out using the program SCHOONSCHIP [15], 
which allows the relevant algebraic manipulations and analytic transformations. The other 

x In Fig. 1, only one of the full set of needed diagrams is presented. 
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important contribution from higher order processes corresponds to the QED rescattering pro- 
cesses in which the scattering pion interacts with the electromagnetic field by exchange of two 
virtual photons. The calculations of the rescattering cross section lead to the analytic results 
presented in [12]. 

These formulae were incorporated into the Fortran code RCFORGV [14]. The RCFORGV 
program calculates more than 100 contributions from the 43 Feynman graphs, that are needed 
to describe the pion scattering by nuclei: 3 diagrams in the Born approximation, one graph 
for the structure effects (polarizability), and 39 diagrams for the radiative corrections. The 
numerical calculations in the RCFORGV program are performed by the Monte Carlo method. 
Double precision numbers are necessary for the calculations of the contribution of the double 
bremsstrahlung process (4). The program package RCFORGV is set up on a IBM-3072 at Tel 
Aviv University. 

The semi-analytical program RCFORGV calculates the radiative corrections of order 0(a) 
to the process (1) for the following cross sections: 

i) for da/du in the w-bins : uj = E\ — E 2 ; 

ii) for d 2 a/dE 2 dQ 2 in the points of (E 2 ,Q 2 ) . 

Here Q 2 = — (pi — p 2 ) 2 , and Ei (E 2 ) are the initial (final) energies of the pions in the 
laboratory system. These cross sections are needed to the analysis of planned experiments. 

In RCFORGV, it is easy to provide a variety of kinematic cuts and criteria for event selec- 
tion, and to take into account the experimental geometry. The following kinematic constraints 
are included in the program: 

1) t = —k 2 = — (pi — p 2 — k') 2 < t for the 4- momentum transfer ; 

2) 9 n < 9 n for the pion scattering angle in the laboratory frame; 

3) uj > uj for the photon detection threshold; 

4) # 77 > # 77 for the angular resolution of the double bremsstrahlung. 

Here t is the maximum value of the squared momentum transfer, fixed by experimental reso- 
lution in t. Below t, Coulomb scattering dominates [7]. Also, 9 n is the maximum of the pion 
scattering angle. The parameters uj and # 77 are important for distinguishing double and single 
bremsstrahlung. If uj < uj or # 77 < 6> 77 , then these two processes are indistinguishable [12]. 

Using RCFORGV, the cross sections at the energies 40 and 600 GeV for the radiative 
scattering of 7r~ mesons on carbon, were calculated in the Born approximation (the Born 
cross section), including polarizability contributions and radiative corrections. Numerical 
calculations were performed by the Monte Carlo method, using 50 000 events for each incident 
energy. The value 40 GeV was chosen to test our new program installation with results given 
previously in [13]. One set of calculations was done for lead, to investigate the Z dependence. 
The cross sections were calculated with kinematic constraints given in Table 1. 

These constraints take into account the geometry and the event selection criteria for ex- 
periments [8] and [10]. Due to kinematics, we chose the parameters of t and uj larger, and 
angles n and # 77 smaller, for energy 600 GeV compared to 40 GeV. 
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Energy 


Parameters 




t 








GeV 


GeV 2 


GeV 


rad 


rad 


40 


2 x 10~ 4 


0.5 


1.5 x 10~ 2 


2.0 x 10~ 3 


600 


4 x 10~ 3 


10.0 


7.0 x 10- 3 


2.0 x 10~ 4 



Table 1: The kinematic constraints used in the program RCFORGV. 

3 The role of radiative correction 

In this section we present the RCFORGV calculations. Following [12], we define the radiative 
correction factor 5(uo) as the ratio of the radiative corrections cross section to the Born cross 
section, as a function of the detected gamma energy u: 

{da/du) rc 

^ = (da/d^om - ( 6 ) 

The results of the RC calculations for the radiative scattering of 40 GeV and 600 GeV n~ 
mesons on carbon nuclei are presented in Fig. 2. One sees that the absolute value of 5(u) is 




Figure 2: Radiative correction 6(u>) for radiative scattering of 40 and 600 GeV it -mesons on 
C 12 nuclei. 

less than 1% over a large energy range. 

The polarizability cross section can be written as the sum of two terms [13]: 

da po1 = -{a n + %)da e + p n da m , (7) 
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where da e and da m give different combinations of electric and magnetic polarizabilities. We 
imposed the condition a n + /3 n = 0. In our analysis, we therefore determine radiative cor- 
rections with respect to one polarizability value = —f3 n . To set the scale, we used the 
values: 



a. 



The ratio of the radiative correction cross section to polarizability cross section is presented 
in Fig. 3. This ratio changes in the interval ~ (—5; 5)% as function of gamma energy uo. The 
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Figure 3: Ratio of the radiative correction to the polarizability cross section for radiative 
scattering of 40 and 600 GeV n~ -mesons on C 12 nuclei. 

formula (5) can be written: 

da tot = da Born + p n da m + da rc . (9) 
To emphasize the fact that radiative corrections can simulate a polarizability, we rewrite (9): 

da tot = da Born + + Pl c )da rn , (10) 

where p r n c = da rc /da m . 

In this analysis, we calculate also the double differential cross section d 2 o-/dE 2 dQ 2 . We 
can transform this cross section to the center mass of the •yn system, as d 2 a / ds\dcos6* . Here, 
s i — (Pi + k) 2 = (p2 + k') 2 and 0* are the 771 energy squared, and the Compton scattering 
angle. At very small —k 2 < t and small 9 n < W we find: 



d 2 a E 1 (s l - mlf{\ - cos6*) d 2 a 



dsidcosO* 



dE 2 dQ 2 



(11) 



5 



The ratio of the polarizability to the Born double differential cross section for different 
scattering angles, as a function of si, is presented in Fig. 4a. In this figure ratios oid 2 a/dE 2 dQ 2 
were taken, but the labeling is according to the CM variables s± and cosO* . It is seen that 
the absolute value of the polarizability contribution increases with energy and can reach a 
value of ~ 30% at backward angles. For a gamma energy u* = (si — m 2 )/2- s /si ~ 150MeV 
and cos9* « —0.5, the contribution of magnetic polarizability to the total cross section of 
reaction (I) is about 7%. At the point s 1 = 0.18GeV 2 (u* « 200MeV) and cos9* = -0.95, 
this contribution attains magnitude of ~ 30%. In the laboratory frame, this kinematical range 
corresponds to the range of detected gamma energies u ~ (0.60 — 0.90)£i, where 



W = ^(i-^)(i-«»n 



(12) 



Here, we point out that data of [8, 9] were taken in this region of uo as well. These 
authors showed that the preferred kinematical range to study the magnetic polarizability in 
the laboratory frame is the range of final gamma energy uo > 0.75i?i. In the center of 
mass frame, this region is defined by the inequalities: si > 0.15GeV 2 and cos9* < —0.75. 
One sees from Fig. 4a that for this kinematical range, the relative contribution of magnetic 
polarizability to the Born cross section is more than 15%. In Fig. 4b, the ratio of the RC to 
polarizability double differential cross section is shown. We see that the relative contribution 
from radiative corrections to the polarizability cross section contribution is less than 5% at 
backward angles for large enough values of si > O.lSGeV 2 . 
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Figure 4: Ratio of a' = d 2 a/dE 2 dQ 2 for: a) the polarizability and the Born cross sections; b) 
the radiative correction and the polarizability cross sections. 



We wish to determine the magnitude of the simulated contribution to the polarizability, 
due to the radiative corrections. We can do this by calculating the average value of the cross 
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section ratio da rc /da m for the gamma energy range to > 0.75E 1 . We have: 

1 n 

= -T, dorc ( u i)/ dom M> ( 13 ) 

where n is the number of u bins. 

For the energies 40 and 600 GeV, we obtained (3 r w c —0.20. The accuracy of our calculation 
may be estimated as: 

where e tot = \Je^ c + e 2 ol . Here, e rc and e po i are the relative uncertainties in the radiative 
correction and polarizability cross section calculations. These uncertainties were calculated 
by the RCFORGV program for each of u bins and we have found that over the gamma energy 
range u> > 0.75E 1 all of e tot are approximately the same . This fact explains the appereance of 
^/n in (14). We have obtained e tot w 0.12, which gives the calculation accuracy 5(/3£ c ) m 0.05. 

In Fig. 5, we show the ratio of the total cross section for reaction (1) to the Born cross 
section. This is given as a function of gamma energy u, for incident pion energies 40 and 600 
GeV. The histogram gives da tot calculated by (9) with /3 n = —5, and the curve corresponds 
to da tot calculated by (10), with j3 n = —5 and (3™ = —0.20. It is clear that both calculations 
are in agreement. The value of (3™ also agrees with the results presented in [13]. The same 
analysis was done assuming that a w — f3 n — 0, and we again find /3£ c ~ —0.20. This last result 
is important, showing that the radiative corrections simulate ff^ ~ —0.20, independent of the 
(3 n value. 




to, GeV 



to, GeV 



Figure 5: The ratio of the total cross section to the Born cross section as a function of detected 
gamma energy u. 

One calculation with (3 n = —5 was carried out to study radiative corrections to pion radia- 
tive scattering on lead. The Born and polarizability cross sections show the Z 2 dependence 
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expected theoretically. The ratio of radiative correction to Born and to polarizability cross 
sections was practically unchanged compared to carbon. This is because 5(u) is practically 
independent of Z. Only the contribution from rescattering has a Z dependence, and this is a 
very small contribution to the total RC. 

We now summarize the main results of our RC study. The radiative corrections simulate 
a polarizability effect at the level of f3™ ~ —0.20. This result is for pion incident energies 
from 40 to 600 GeV, and for nuclei from carbon to lead. We find for (3™ ~ _ .20 ± 0.05 at 40 
and 600 GeV, for /3 n from the interval (—7; —2), the relative contribution to the polarizability 
from radiative corrections equals (3^/ (3^ ~ (3 — 10)%. Therefore, we give a simple lowest order 
method of obtaining the polarizability from experimental data. We write the cross section as 
a sum of Born and polarizability cross sections: 

da exv = da Born + fi exp do m . (15) 

We now find the best fit /3% xp for da exp using RCFORGV, or any other convenient fit pro- 
gram. This experimental /3%. xp must be corrected by [3^°. We have by comparison to (10), the 
relationship: j3* xp = (3 n — (3™. Thus a value fi exp = —5.20 would correspond a polarizability 
fin = —5. A more exact determination would require analysis of the data with RCFORGV. 



4 Different Born cross section calculations 

In this section, we briefly describe two methods for the Born cross section calculation. The first 
is an exact calculation performed with the help of the RCFORGV program, and the second 
is the equivalent photon approximation (EPA) (Weizsacker- Williams's method) [16, 17]. 

The exact calculation method consists of an accurate counting of the contributions from 
the all QED diagrams in the lowest order of a. This corresponds to the scattering of a spinless 
pointlike particle by a Coulomb center with hard photon emission. The exact calculation takes 
into account the magnitude of the photon virtuality k 2 . The exact expression of the Born cross 
section of the process (1) has the following form [12]: 

d 2 a Born tfrf - rff 



dE 2 dQ 2 f 2 J tmin t 2 ' S ° {t) - (16) 

The explicit form of the function S°(t) and the kinematical limits t m i n ^ max are given 
in [11, 12]. Here we emphasize that the exact Born cross section calculation includes the 
contribution of the logarithmic term ~ ln(t max /t min ) (leading logarithmic approximation), 
and the additional contribution of the next term which is proportional to ~ (l/t min — l/t max ). 

The EPA method is based on the assumption that the photon virtuality is very small, 
—k 2 < t << m 2 . Therefore, one ignores the small deviation from the real photon, and one 
takes only the leading logarithmic term. In this case, the cross section of the reaction (1) can 
be expressed via the cross section for the 77? scattering (2): 

da vz ^*nz = n(si)dsid<r 7W (si). (17) 

Here, d<T 77r (si) is the cross section of the Compton scattering of the real photon, n(si) is the 
equivalent photon density . 

In the EPA, the cross section for the reaction (1) in the center of mass frame can be written: 



d (J Ot Z , tm.n.cr. tm.ir) ^ d(T> 



ds\dcos6* it si — m 2 



1 -max , v mm 

In 1 



777 



dcosQ 



(18) 
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where d<7 77r (si) / dcosO* is the unpolarized differential cross section of the reaction (2). 
The Born cross section of the 77r Compton scattering has the form [18] : 



da, 



Born 

77T 



where ti = (pi 



dcosO* 
P2) 2 ,u l = (p 



l + [-l + 



(si - ml)(u! - ml)- 



(19) 



k')' 
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Sj, squared of ( y+Jl ) energy , ( GeV 2 ) 



+ cos*0 = -.95 
□ cos'O = -.55 
A COS*0 = -.15 



Figure 6: Comparison of two Born cross section calculations: 
method (a' = d 2 a/dE 2 dQ 2 ). 



exact calculation and EPA 



A comparison between the two Born cross section calculations was carried out in order 
to estimate their applicability. In Fig. 6a the difference between the two double differential 
cross sections d 2 a/dE 2 dQ 2 , is presented. The superscript Born refers to the exact Born cross 
section performed by RCFORGV program, and WW to the Born cross section calculated by 
the EPA method. As in Fig. 4, ratios of d 2 a/dE2dQ 2 were taken, but the labeling is according 
to the CM variables si and cos9*. The cross section difference is given with respect to the 
exact Born cross section. 

From Fig. 6a, we can see that the EPA method agrees with the exact Born calculation 
with a precision better then 3 % in the kinematical region si > 0.13GeV 2 and at backward 
scattering angles cosO* < —0.5. From Fig. 6b, the cross sections difference divided by the 
polarizability cross section is displayed. In the kinematical range limited by the inequalities 
si > Q.lSGeV 2 and cosO* < —0.5, both Born cross section calculation methods are in agree- 
ment and can be applied equally well. This kinematical region was also preferred for the study 
of the magnetic polarizability (see section 3). 

The validity of the EPA method was discussed in recent papers (see for example [19]- [21]) 
shown that the EPA for the other processes can lead to large errors, if used outside of the 
region of applicability and must be used judiciously or improved in order to get correct results. 
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5 Conclusions 



We used the semi- analytical program RCFORGV to evaluate the contribution of the radiative 
corrections to the total cross section of the pion scattering by a Coulomb center with hard 
photon emission. We showed that the radiative corrections can simulate the polarizability 
effect and the average effect (3'^ c = (—0.20 ± 0.05) x 10~ 43 cm 3 was obtained. We showed 
that the preferred kinematical region to investigate the magnetic polarizability cross section 
in the laboratory frame is the range of final gamma energy uu > 0.75-Ei . The corresponding 
range in the center of mass frame is defined by the following inequalities : si > 0.15GeV 2 
and cosO* < —0.75, where s\ and 9* are the energy squared and the polar scattering angle. 
The cross sections of the reaction (1) computed in the Born approximation using the exact 
calculation method and the equivalent photon approximation are in the agreement for hard 
emitted photons u > 0.75£i . In the range of the photons uo < 0.6-Ei, we must use the exact 
calculation method. The suggested method allows to take into account the radiative correction 
effects with high accuracy. 
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